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We present searches for rare or forbidden charm decays of the form X^ -^ h (.^ t^'^ , where 
X^ is a charm hadron (-D^, Dt , or At), h is a pion, kaon, or proton, and t^' is an electron or 
muon. The analysis is based on 384 fb~^ of e+e~ annihilation data collected at or close to the T{AS) 
resonance with the BABAR detector at the SLAC National Accelerator Laboratory. No significant 
signal is observed for any of the 35 decay modes that are investigated. We establish 90% confidence- 
level upper limits on the branching fractions between 1 x 10~^ and 44 x 10~^ depending on the 
channel. In most cases, these results represent either the first limits or significant improvements on 
existing limits for the decay modes studied. 

PACS numbers: 11.30.Fs,11.30.Hv,13.30.Ce,13.20.Fc 



I. INTRODUCTION 



We present searches for charm hadron decays that are 
either forbidden or heavily suppressed in the Standard 
Model (SM) of particle physics. The decays are of the 
form ^ A+ -^ h'^i^&^^ , where X^ is a charm hadron 
(D+, I?+, or A't), and i'^'^^ is an electron or muon. For 
£*+ and 73+ modes, h^ can be a pion or kaon, while 
for A't modes it is a proton. Decay modes with oppo- 
sitely charged leptons of the same lepton flavor are exam- 
ples of flavor-changing neutral current (FCNC) processes, 
which are expected to be very rare because they cannot 
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occur at tree level in the SM. Decay modes with two 
oppositely charged leptons of different flavor correspond 
to lepton-flavor violating (LFV) decays and are essen- 
tially forbidden in the SM because they can occur only 
through lepton mixing. Decay modes with two leptons 
of the same charge are lepton-number violating (LNV) 
decays and are forbidden in the SM. Hence, decays of 
the form A+ — )■ h^H.'^i'^'^^ provide sensitive tools to in- 
vestigate physics beyond the SM. The most stringent ex- 
isting upper limits jli-SJ on the branching fractions for 
X+ -^ /i±£=F£(')+ decays range from 1 to 700 x 10"^ and 
do not exist for most of the yl+ decays. 

FCNC processes have been studied extensively for K 
and B mesons, in K^ — K^ and B^ — B^ mixing, and 
in rare FCNC decays such as s ^' d£^i^, b -^ sj, and 
b — !■ s£^£^ decays. The results agree with expectations 
within the framework of the SM Q . There are ongoing 
efforts to improve the measurements and the theoreti- 
cal predictions, and to measure new effects, such as CP 
violation, in FCNC processes. 

The recent observation of D'^ — D*' mixing [3-[3 has in- 
creased interest in FCNC processes in the charm sector. 
Of particular interest is the source of D^ — D^ mixing. If 





FIG. 1: Standard model short-distance contributions to the 
c -> ut^(.~ transition. 



II. OVERVIEW 

The searches use data collected with the BaBAR de- 
tector as described in Section Hill Candidate decays are 
formed from two tracks identified as leptons (£^^'^) and 
one track identified as a hadron {h). Background events 
arise primarily from semileptonic B and charm decays 
and from radiative QED events with converted photons. 
After applying basic kinematic selection criteria, we dis- 
criminate further between signal and background with 
a decay-mode-dependent likelihood ratio calculated from 
the measured momentum and flight length of the charm 
hadron candidate and the total energy detected in the 
event. The invariant mass distributions of the selected 
/i£^('' candidates are fit to extract the number of signal 
events, the number of combinatorial background events, 
and the number of background events due to nonleptonic 
charm decays with hadrons in the final state misiden- 
tified as leptons. Upper limits are determined for the 
decay rate to each final state with respect to the known 
decay rate for a hadronic three-body decay with similar 
kinematics so that many systematic uncertainties cancel 
in the ratio. 



the mixing is due to physics beyond the SM, it could also 
give rise to measurable effects in FCNC charm decays. In 
the SM, the FCNC decays X+ -^ h+£+£- are expected 
to be heavily suppressed due to cancellations of am- 
plitudes thr oug h the Glashow-Iliopoulos-Maiani (GIM) 
mechanism [lOJ. For example, the c -^ ut^t~ transi- 
tions illustrated in Fig. [T] yield branching fractions for 
D ^ X„£+£- of 0(10-8) [HI, [13. These decays are 
masked by the presence of long-distance contributions 
from intermediate vector resonances such as 13 — > X^y^ 
V -^ £'^£~, which are predicted to have branching frac- 
tions of 0(10-^) [111, [l2|- The effect of these resonances 
can be separated from those due to short-range processes 
by applying selection criteria on the invariant mass of the 
£'^£~ pair. In radiative charm decays, c -^ 1*7, uncertain- 
ties in calculating the long-distance terms make it impos- 
sible to study the underlying short-distance physics [13|. 

The impact of several extensions of the SM on D -^ 
Xu£^£^ decay rates has been estimated [lil, [ij, [i3 J the 
largest effect not already ruled out is expected in certain 
R-parity violating supersymmetric models. Depending 
on the size of the R-parity violating couplings, branching 
fractions of up to 0(10-^) for different D -^ Xu£+£~ 
decays are possible. 

The only long-distance amplitudes relevant at the cur- 
rent experimental sensitivity are from £)+ — > 7r+0 and 



Df 



-^ n^ 



^(K) 



decays. The product of branching fractions, 



n+cj))B{(j) -^ £+£-), is f* 2 X 10"^ for D+ and 

sa 1 X 10^^ for Df 16]. In this analysis, we measure the 
total decay rate of Dt. -^ tt~^£'^£~ excluding the £+£^ 
mass region around the (p resonance. 



III. THE BABAR DETECTOR AND DATASET 

The BaBAR detector was operated at the PEP-II 
asymmetric-energy storage rings at the SLAC National 
Accelerator Laboratory. The data sample comprises an 
integrated luminosity of 347 fb" collected from e+e~ 
collisions at the T(4S') resonance and 37 fb" collected 
40 MeV below the r(4S') resonance. 

The BaBAR detector is described in detail else- 
where [ITI . The momenta of charged particles are mea- 
sured with a combination of a five-layer silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH), both 
in a 1.5 T magnetic field produced by a solenoid. The 
resolution of the transverse momentum (px) is measured 
to be (j{pt)/pt = 0.0013(pt/ GeV/c) ® 0.0045. A detec- 
tor of internally reflected Cherenkov radiation (DIRC) 
is used for charged particle identiflcation. Pions, kaons, 
and protons are identified with likelihood ratios calcu- 
lated from dE/dx measurements in the SVT and DCH, 
and from the observed pattern of Cherenkov light in the 
DIRC. For the selection criteria used in this analysis, 
hadron identification efficiencies are approximately 98%, 
87%, and 82% for pions, kaons, and protons, respectively. 
A finely segmented CsI(Tl) electromagnetic calorimeter 
(EMC) is used to detect and measure the energies of pho- 
tons and neutral hadrons, and to identify electrons. For 
electrons, energy lost due bremsstrahlung is recovered 
from deposits in the EMC. 

The EMC information, Cherenkov angle, and dE/dx 
measurements are combined to define a likelihood ratio 
used to select electrons with a selection efficiency above 
90% for electrons with a laboratory momentum above 
0.5 GeV/c and with less than 0.1% probability to misiden- 



tify a hadron as an electron. The instrumented flux re- 
turn (IFR) contains resistive plate chambers and limited 
streamer tubes for muon and long-lived neutral-hadron 
identification. Variables characterizing track measure- 
ments in the IFR and the energy deposition in the EMC 
are combined in a neural network to select muon can- 
didates. The muon identification efficiency is about 60% 
for muons with a laboratory momentum above 1.5 GeV/c, 
but decreases rapidly for lower momenta. The probabil- 
ity to misidentify a pion as a muon is about 1.5% for 
most of the relevant momentum range. 

Event simulations are performed using the EvtGen [l8| 
Monte Carlo (MC) generator with a full detector simula- 
tion based on GEANT4 ^. Signal and A+ -^ pR- Tr+ 
MC events are generated with a three-body phase-space 
distribution, while Dt-. -> tt^<J) MC events are gener- 
ated with a Breit-Wigner resonance shape for the (p decay 
and an angular distribution appropriate for the P-wave 
(/> decay. All signal events are simulated as cc contin- 
uum events. The distributions of the magnitudes of the 
charm-hadron momenta measured in the e~^e~ center-of- 
mass frame (p* ) are found to differ between simulated cc 
continuum events and data. In particular, the mean p* 
for D+ mesons is lower in MC, while the mean p* for 
/!+ baryons is higher. To correct for this, simulated cc 
events are weighted to yield the same p* distributions as 
those observed for large samples of reconstructed charm 
decays in data. The event weights lie between 0.6 and 1.8. 
Simulated events are also weighted to match the particle 
identification probabilities measured in control samples 
in data. Samples of simulated cc and qq (q — u, d, s) 
continuum events and BB events, corresponding to 1.4 
to 5 times the recorded data sample, are used to study 
background contributions. 



IV. ANALYSIS PROCEDURE 
A. Initial Signal Selection 

Charm hadron candidates are formed from one track 
identified as either a pion, kaon, or proton (h) and two 
tracks, each of which is identified as an electron or a 
muon {££'-''> ) . The total charge of the three tracks is re- 
quired to be ±1. For three-track combinations with a 
pion or kaon track, the M£^'' invariant mass is required 
to lie between 1.7 and 2.1 GeV/c^; for combinations with 
a proton, the invariant mass is required to lie between 
2.2 and2.4GeV/c2. 

The combinatorial background at low p* is very large 
and we therefore select charm hadron candidates with 
p* greater than 2.5 GeV/c. Because the p* of charm 
hadrons produced in B decays is kinematically limited 
to be less than about 2.2 GeV/c, the signal candidates 
with p* greater than 2.5 GeV/c are dominated by hadrons 
from continuum production. The main backgrounds re- 
maining after this selection are QED events and semilep- 
tonic B and charm decays, particularly events with two 



semileptonic decays. 

The QED events are mainly radiative Bhabha, initial- 
state radiation, and two-photon events, which are all 
rich in electrons. These events are easily identified by 
their low multiplicity and/or highly jet-like topology. We 
strongly suppress this background by requiring at least 
five tracks in the event and that the hadron candidate be 
inconsistent with the electron hypothesis. 

We suppress the background from semileptonic B and 
charm decays by requiring the two leptons to be consis- 
tent with a common origin. This is achieved by requiring 
that the probability of the x^ for the vertex fit be greater 
than 0.001 and that the distance of closest approach be- 
tween the two lepton candidates be less than 250 /xm. 

For low e+e" invariant mass there is a significant back- 
ground contribution from photon conversions and tt'^ de- 
cays to e^e^7. These are both removed by requiring 
m(e+e-) > 200MeV/c2. 

For the Dfs -^ 7r+£+f^ decay modes, we exclude 

events with 0.95 < m{e+e-) < 1.05GeV/c2 and 0.99 < 
m{^^ p.~) < 1.05GeV/c^ to reject decays through the 
resonance. The excluded regions for the two decay modes 
differ due to the larger radiative tail in the 7r+e^e~ de- 
cay mode. In order to perform cross-checks, we also se- 



lect candidates for D 



(s) 



—^ 7r0, (/) — i> £~^£ decays by 

requiring 0.995 GeV/c^ < m{e+e-) < 1.030 GeV/c^ or 
1.005 GeV/c2 < m{fi+fi-) < 1.030 GeV/c^. 



B. Likelihood Selection 

After the initial event selection, significant combinato- 
rial background contributions remain from semileptonic 
B decays and other sources. These background sources 
are studied with samples of candidates in MC and in side- 
bands of the h££^'' invariant mass distribution in data. 
The sidebands are defined to be 0.1 GeV/c^ wide regions 
just below and above the signal regions defined above. 
The final candidate selection is performed by forming a 
likelihood ratio Re and requiring the ratio to be greater 
than a minimum value i?™'". The likelihood ratio is de- 
fined as 



Rc{x) = 



Il^^^sA^^) 



a^s 



U^'Pb4^^) 



(1) 



where Xi is the ith. discriminating variable, and Vs,i and 
Vb,i are the signal and background probability density 
functions (PDF) for the variable Xi. Correlations be- 
tween variables are ignored in the likelihood ratio as they 
are found to be small. The likelihood ratio peaks near 1 
for signal and near for background. The signal PDFs 
are obtained by fitting the distribution of Xi for simulated 
signal events, while the background PDFs are obtained 
from fits to the distribution of Xi for candidates in the 
h££^'^ invariant-mass sidebands in data. The PDFs are 
defined by combinations of polynomial, Gaussian, and 
exponential functions found empirically to provide good 



descriptions of the signal MC and data sidebands. The 
parameter values in the PDFs are common across signal 
modes in which data or MC studies show the distribu- 
tions to be consistent. 

The following three discriminating variables are used 
in the likelihood ratio. 

• Charm hadron candidate p* . 

The calculated p* for hii'^'^ candidates in which a 
lepton candidate is from a semileptonic B decay 
can be larger than 2.2 GeV/c, but the distribution 
falls rapidly with increasing p* . 

• Total reconstructed energy in the event. 

Since neutrinos from semileptonic B decays are not 
directly detected, the total observable energy in 
these events is less than the sum of the beam ener- 
gies (12GeV). We calculate the total energy for 
each event as the sum of the energies of all re- 
constructed tracks (assuming each track to be a 
charged pion) and neutral EMC clusters. 

• Flight length significance. 

The flight length significance is the ratio of the 
signed flight length to its uncertainty, where the 
signed flight length is the scalar product between 
the direction of the charm-hadron candidate and 
the position vector that points from the beam spot 
to the charm-hadron decay vertex. This variable 
has the most discrimination power for the long- 
lived D+ meson and the least for the shorter-lived 
yl+. This variable is effective for suppressing non-_B 
combinatorial background. 

Distributions of the three discriminating variables and 
of the likelihood ratio are shown in Fig. [2] for _D+ -^ 
TT+e+e^ candidates in hi£^'' invariant-mass sidebands 
(background) and in signal MC. 

The minimum likelihood ratio value _R™'" is chosen in- 
dependently for each mode to provide the lowest expected 
upper limit on the branching ratio, as calculated from 
the simulated signal efficiency and the expected number 
of background events; the latter is estimated from the 
Mi^'' invariant-mass sidebands in data. In cases in which 
more than one candidate from the same event passes all 
selection criteria, including the likelihood-ratio require- 
ment, the candidate with the highest p* is retained. The 
final signal selection efficiency lies between 0.5% and 7%, 
depending on the signal mode. 



C. Fit Procedure 

Extended, unbinned, maximum-likelihood fits are 
applied to the invariant-mass distributions for the 
h^ i"^ £'-''> ~^ candidates. The PDF we use for signal events 
is the so-called Crystal Ball function [20], which has an 
asymmetric component to describe the radiative tail in 
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FIG. 2: Distributions of the three discriminating variables 
and the likelihood ratio, for D'*' -+ n'^e^e^ candidates in sig- 
nal MC (black solid curve) and in the h£r'' invariant-mass 
sidebands from data (red dashed curve): a) Center-of-mass 
momentum p* of the charm-hadron candidate; b) Total en- 
ergy in the event; c) Flight length significance; d) Likelihood 
ratio (defined in Eq. [l]) calculated with the three discriminat- 
ing variables. The signal and background distributions are 
normalized to the same area. 



the mass distribution: 



PcB{m;fi,a,a,n) 

p 2 V CT -* 



na — a^a — a{m—fi) 



a m > fi — aa, 
if TO < /x — aa. 



The four parameters /i, a, a and n are determined from 
fits to signal MC candidates and are fixed to these values 
in the fits to data, with only the overall normalization as 
a free parameter. The width of the Gaussian component 
(ct) is found to lie between 5 and lOMeV/c^, depending 
on the decay mode. 

The invariant mass distributions of the combinatoric 
background events for the signal modes are described by 
first-order polynomials. The background slope is a free 
parameter in all fits. 

An additional background arises from nonleptonic 
charm decays in which two hadrons are misidentified as 
leptons. This background component is almost negligi- 
ble in the signal modes with electrons and is therefore 
only included in decay modes with two muons and for 






K- 



+^+ 



/i ' e 



where there is a large D+ — > K^tt^tt^ 
background. The shape of this background is obtained 
from MC samples of major hadronic three-body charm 
decays. Each event is weighted according to the proba- 
bility of misidentifying a pion as a lepton. The misiden- 
tification probability is measured from data as a function 
of momentum and angle with samples of D" — > K^tt'^ 
candidates. The misidentified nonleptonic charm decays 
are reconstructed at slightly lower /i+/z"'"/j,~ {K~ n~^e~^) 
mass than the signal events. The peak mass is shifted 
by about I5MeV/c^, which is sufficient separation from 



TABLE I: Branching fractions for the charm decays used for 
normalization llql. 



Decay mode 



Branching Fraction 



At 






(2.72±0.13) X 10"^ 

(2.32±0.14) X 10"^ 

(5.0±1.3) X 10'^ 



TABLE II: Fitted signal yields for the normalization modes 
and signal efficiencies estimated from MC simulations. Only 
statistical uncertainties are quoted. 



Decay mode 



At 



> -K (pKK 



N,i 



106 800 ± 500 

338 900 ± 900 

488 700 ± 2 100 



Efficiency 



(15.44 ± 0.07)% 
(15.29 ±0.07)% 
(11.99 ±0.04)% 



the signal peak for this background yield to be deter- 
mined by the likelihood fit to data, without reliance on 
MC predictions of this yield. 



V. BRANCHING RATIO NORMALIZATION 

The measured signal yields are converted into branch- 
ing ratios by normalizing them to the yields of known 
charm decays. We choose normalization modes with 
kinematics similar to the kinematics of the signal decays 
so that most systematic effects not related to particle 
identification cancel in the branching ratio. For the D^ 
and Z?+ mesons, we use decays to 7r"*"0 as normalization 
modes, where the cj) decays to K^K~ . For the yl+, we 
use At -^ pA'^7r+ as the normalization mode. The mea- 
sured branching fractions for these modes are listed in 
Table HI We use the abbreviation Df^-^ — >■ t:^4>kk to 

denote D 



(s) 



K^K~ decays. 



A. Event Selection 

The same selection criteria are applied for the normal- 
ization modes as the signal modes, except for the lep- 
ton identification and likelihood ratio requirements. In- 
stead, each daughter candidate in the Dt^ -^ t^^(J}kk 

and At — > pK^ir^ decay modes is required to be consis- 
tent with the kaon, pion, or proton hypothesis, as appro- 
priate. For the Dt. -^ tt~^4'kk decay modes, we further 
require the invariant mass of the kaon pair to lie within 
15MeV/c^ of the world- average mass [la |. 



The D 



(s) 



'^^4'KK samples include a small compo- 



nent of non-resonant or S-wave Dt-^ decays, while the 
branching fractions in Table U are extracted from de- 
cay amplitude analyses of Dalitz plot distributions for 



D 



is) 



TT^K^K decays and therefore correspond only 
to the resonant component. This component is estimated 
in our MC and data samples by projecting out the P-wave 
component by weighting each event with a factor that in- 
cludes the reciprocal of its reconstruction efficiency and 
a normalized L=2 Legendre polynomial function of the 
K~^K~ helicity angle. The fractions of D'^ and Df de- 
cays found to proceed through a P-wave are 94% and 
93%, respectively, and are used to correct the fit yields. 




B. Fit Results 



2.28 2.29 2.3 2.31 2.32 

pK*jt" invariant mass [GeV/c^] 



FIG. 3: Invariant mass distribution for a) n^fjiK+K- ^nd b) 
pKir candidates. The solid lines are the results of fits to 
double- Gaussian distributions for signals and a second-order 
polynomial for the background (dashed line). 



For the normalization modes, radiative effects are neg- 
ligible and we use the sum of two Caussian distributions 
with a common mean to describe each of the 13+ , Df , 
and At signal invariant-mass distributions. All signal pa- 
rameters are free in the fits to the normalization modes 
in data. The combinatorial background is described by 
a second-order polynomial. The invariant-mass distribu- 
tions for the normalization decay modes and the corre- 
sponding fit results are shown in Fig. [31 The fitted signal 
yields are listed in Table HIl where we also list the effi- 
ciencies estimated from signal MC. 



VI. SYSTEMATIC UNCERTAINTIES 

Most systematic effects are expected to cancel in the 
branching ratio since they affect both the signal and nor- 
malization modes. Therefore, the main systematic un- 
certainties arise from differences in selection, acceptance, 
and decay kinematics. Table IIIII gives a summary of all 
systematic uncertainties related to the branching ratio 
calculation. An additional systematic uncertainty is as- 
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TABLE III: Summary of the multiplicative systematic un- 
certainties on the branching ratio for each signal decay mode 
(listed in the first column). In the second column, we give the 
fractional uncertainty due to the normalization mode. The 
next three columns list the fractional uncertainty due to the 
MC statistical error, uncertainties in particle identification 
efficiencies (PID), and in the likelihood-ratio efficiency (LR). 
The total systematic uncertainty, given in the final column, 
is the sum of the individual errors in quadrature. 





Norm. 










Decay mode 


Mode 


MC Stat. 


PID 


LR 


Total 


D+ -^n+e+e- 


2.6% 


1.8% 


2.3% 


9.8% 


10.6% 


D+ -^ TT + fl+fl- 


2.6% 


4.6% 


8.3% 


6.1% 


11.6% 


D+ -^ -K+e+^i- 


2.6% 


3.1% 


4.4% 


9.8% 


11.5% 


D+ -^ TT+fi+e- 


2.6% 


2.2% 


4.4% 


7.6% 


9.4% 


D+ -^ TT+e+e- 


2.1% 


1.1% 


2.3% 


0.9% 


3.4% 


D+ -^ TT+n+fi- 


2.1% 


2.7% 


8.3% 


0.9% 


9.0% 


D+ -^ TT+e+fj.- 


2.1% 


2.6% 


4.4% 


2.9% 


6.2% 


D+ -^ 7r+/i+e- 


2.1% 


4.0% 


4.4% 


7.2% 


9.6% 


D+ -^ K+e+e- 


2.6% 


1.4% 


2.8% 


5.5% 


6.8% 


D+ -^ K+n+fi- 


2.6% 


6.4% 


8.5% 


4.4% 


11.8% 


D+ -^ K+e+n' 


2.6% 


2.9% 


4.7% 


5.8% 


8.4% 


D+ ^ K+ti+e- 


2.6% 


3.1% 


4.7% 


5.1% 


8.0% 


D+ -^ K+e+e- 


2.1% 


1.5% 


2.8% 


2.0% 


4.3% 


D+ -> K+n+fi- 


2.1% 


3.3% 


8.5% 


0.9% 


9.4% 


Dt -^ K+e+^i- 


2.1% 


2.2% 


4.7% 


2.0% 


5.9% 


Df -^ K+n+e- 


2.1% 


2.0% 


4.7% 


1.6% 


5.7% 


At -> pe+e~ 


3.4% 


1.4% 


2.0% 


5.4% 


6.8% 


At -> p^J,^^J,~ 


3.4% 


3.2% 


8.2% 


3.4% 


10.0% 


At -> pe+^~ 


3.4% 


3.4% 


4.3% 


5.6% 


8.5% 


At -> pfi^e~ 


3.4% 


2.8% 


4.3% 


5.1% 


8.0% 


D+ -^n-e+e+ 


2.6% 


1.3% 


2.3% 


5.7% 


6.8% 


D+ -^ TT~fl + ^J,+ 


2.6% 


3.5% 


8.3% 


5.1% 


10.7% 


D+ -^TT'fi+e+ 


2.6% 


1.6% 


4.4% 


4.6% 


7.1% 


Dt -^ 7r-e+e+ 


2.1% 


0.8% 


2.3% 


1.2% 


3.4% 


Dt -^ 7r^/i+/i+ 


2.1% 


2.2% 


8.3% 


1.2% 


8.9% 


D+ -^ TT~fj.+ e+ 


2.1% 


2.4% 


4.4% 


1.2% 


5.6% 


D+ ^K-e+e+ 


2.6% 


1.4% 


2.8% 


5.7% 


7.0% 


D+ -> K-fi+ij.+ 


2.6% 


3.2% 


8.5% 


2.7% 


9.8% 


D+ -^ K-fi+e+ 


2.6% 


2.4% 


4.7% 


3.9% 


7.1% 


D+ -^ K-e+e+ 


2.1% 


1.0% 


2.8% 


1.4% 


3.9% 


D+ -i- K-ij,+tj.+ 


2.1% 


2.8% 


8.5% 


0.6% 


9.2% 


D+ -^ K-fi+e+ 


2.1% 


1.9% 


4.7% 


0.6% 


5.5% 


At -^ pe+e+ 


3.4% 


1.0% 


2.0% 


5.4% 


6.8% 


At -> m'^fJ-'^ 


3.4% 


2.6% 


8.2% 


1.7% 


9.4% 


At -^ pfi+e+ 


3.4% 


1.6% 


4.3% 


3.5% 


6.7% 



sociated with the estimate of the signal yield. 

Systematic uncertainties related to the signal PDF pa- 
rameters are investigated in two ways. First, the PDF 
parameters for data and MC are compared in the fits 
to the normalization modes. Differences can be due ei- 
ther to general data-MC tracking differences or to dif- 
ferences between the simulated charm hadron mass and 
the actual mass. Second, fits to the invariant mass of 
J/V' — >■ i'^i~ candidates from inclusive B decays are com- 
pared between data and MC. The second comparison is 



sensitive to effects associated with lepton reconstruction 
and radiative tails. Based on these studies, the mean sig- 
nal mass is varied by up to 2.5 MeV/c^, depending on the 
decay mode. Based on the same studies, the widths of 
the signal PDFs are varied by ±5%. For the background 
shape assumption, the signal fits are repeated using a 
second-order polynomial as the background PDF instead 
of the nominal first-order polynomial. The PDF varia- 
tion giving the worst upper limit on the branching ratio 
is used as an estimate of the systematic uncertainty. 

For the normalization modes, the statistical uncertain- 
ties from the fits, the MC statistical errors, and uncer- 
tainties from the signal and background shapes are all 
about 1% or less. The main systematic uncertainty is 
due to intermediate resonances in the three-body decays. 
For the Dt. modes, where we select K'^K' pairs con- 
sistent with the (f) mass, this uncertainty is estimated by 
varying the K~^K~ mass interval in order to alter the pu- 
rity of the accepted (j> candidates. For the At mode, the 
uncertainty is estimated from variations in the efhciency 
as a function of the pK~ and pn'^ invariant mass. For 
all three normalization modes the estimated uncertainty 
is approximately 2%. 

The particle identification efficiencies have associated 
systematic uncertainties of 0.5% for each pion, 0.7% for 
each kaon, 0.9% for each electron, and 4% for each muon. 
We do not evaluate a systematic uncertainty for the pro- 
tons since, for the Ac decays, both the signal and nor- 
malization modes contain a proton and the uncertainty 
therefore cancels. Uncertainties for particles of the same 
type are added linearly, while those for different types of 
particles are added in quadrature. 

The systematic uncertainty on the efficiency of the 
likelihood-ratio requirement is estimated by applying the 
same likelihood ratio selection to events in the normal- 
ization mode. The efficiency of this selection for the 
normalization-mode decays is not expected to be the 
same as for signal modes due to different kinematics. 
However, for the normalization modes, the efficiency can 
be measured for both data and for MC simulation; the 
difference is used as the systematic uncertainty. The 
largest variations are found for decay modes with the 
most stringent likelihood-ratio requirements. 

In the calculation of the signal efficiency we assume 
that the decays populate three-body phase-space uni- 
formly. The selection efficiency has some dependence on 
where the decay lies in the Dalitz plot. Ignoring the re- 
gions we explicitly remove in the selection, the efficiency 
varies by less than 25% around the average as a function 
of the i^i'^'''^ invariant mass, with the lowest efficiency 
at low dilepton mass. This model dependence is not in- 
cluded in the systematic uncertainty. 



VII. RESULTS 

The hi£^'' invariant-mass distributions for signal can- 
didates in all 35 decay modes are shown in Figs. HHHl 
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The signal yields obtained from the unbinned likelihood 
fits are listed in Table IIVI with statistical and system- 
atic uncertainties. Only systematic uncertainties asso- 
ciated with the signal and background PDFs are in- 
cluded in the systematic uncertainty for the yields. The 
curves representing the fits are overlaid in the figures. 
The most significant signal is seen in the distribution 
for Af -^ p^~^ii~; the signal yield has a statistical-only 
significance of 2.6a as determined from the change in 
log-likelihood with respect to zero assumed signal events. 
With 35 different measurements, a 2.6a deviation is ex- 
pected with about 25% probability. 

We calculate upper limits on the branching ratios 



is) 



B{D+ ^ K^i^li^''^+)/B{D+ -^ ^+0), and 



(s) 



is) 



B{At ^ ^p^tf(.^''^+)/B{At ^pK-Ti+) 

at 90% confidence level (CL). The upper limits are set 
using a Bayesian approach with a flat prior for the event 
yield in the physical region. The upper limit on the sig- 
nal yield is defined as the number of signal events for 
which the integral of the likelihood from zero events to 
that number of events is 90% of the integral from zero to 
infinity. The systematic uncertainties are included in the 
likelihood as additional nuisance parameters. For com- 
parison with previous measurements, the upper limits on 
the branching fractions at 90% CL, calculated using the 
data of Table HI are also given. 

For 32 of the 35 decay modes, this analysis is more 
sensitive than existing measurements. In most cases, the 
improvement is significant (factor of 2 to 60). The largest 
improvements are seen for the LFV decays, which are 
all improved by at least a factor of 10. The only A'^ 
decay with a pre-existing limit is A'^ — > pfi^ fi~ , which we 
improve by roughly a factor of 8. For all other Af decays 
this paper presents the first limits. The only modes that 
do not provide a more sensitive limit are D^ — ^ Tr^e+e^, 
D~^ — >■ Tr~^ fi'^ n~ , and Df — >■ n^ii'^fi^ where existing 
limits [l|-yl are about a factor of two lower than those 
presented here. 



As a cross check, we also examine D 



D 



r+. 



is) 



and 



, > -^ n v^^+^- events with dilepton invariant masses 
in the (j) signal region defined in Sect. IIV Al The in- 
variant mass distributions are shown in Fig. 1101 Signals 
with a statistical significance of at least 3a are seen for 
all decays except for Z3+ -^ 7r+0^+^-. The selection 
for _D+ and Df candidates differ through the likelihood- 
ratio criteria, which are optimized separately; however, 
signals for both hadrons can be seen with either selection. 
The mass distributions are therefore fit allowing for both 
a D~^ and D^ signal, but only the signal yield for the 
hadron for which the likelihood ratio was constructed is 
used. Table IVl gives the fit yields. It also shows the ex- 
pected yield assuming 6(0 -^ £'^£~) = B{(f> -> e+e^) — 
(2.95 ±0.03) X 10"'' [ij]. The fit yields are in good agree- 
ment with expectations. 



TABLE IV: Signal yields for the fits to the 35 X+ -5> 
h (.^ t-''^ event samples. The first error is the statistical un- 
certainty and the second is the systematic uncertainty. The 
third column lists the estimated signal efficiency. The fourth 
column gives for each signal mode the 90% CL upper limit 
(UL) on the ratio of the branching fraction of the signal mode 
to that of the normalization mode (BR). The last column 
shows the 90% CL upper limit on the branching fraction for 
each signal mode (BF). The upper limits include all system- 
atic uncertainties. 









BRUL 


BF UL 




Yield 


Eff. 


90% CL 


90% CL 


Decay mode 


(events) 


(%) 


(10-*) 


(10-') 


D+ ^TT+e+e" 


-3.9 ±1.6 ±1.7 


1.56 


3.9 


1.1 


D+ -> 7r+/i+^" 


-0.2 ±2.8 ±0.9 


0.46 


24 


6.5 


D+ ^ Tv+e+fj.- 


-2.9 ±3.4 ±2.4 


1.21 


11 


2.9 


D+ -^ 7r+/i+e- 


3.6 ±4.3 ±1.3 


1.54 


13 


3.6 


Dt ^TT+e+e" 


8 ± 34 ± 8 


6.36 


5.4 


13 


Dt -> 7r+/i+^" 


20 ± 15 ± 4 


1.21 


18 


43 


Dt -^ Tv+e+fj.- 


-3 ±11 ±3 


2.16 


4.9 


12 


Dt -> 7r+/^+e- 


9.3 ±7.3 ±2.8 


1.50 


8.4 


20 


D+ -^ K+e+e- 


-3.7 ±2.9 ±3.3 


2.88 


3.7 


1.0 


D+ -> K+fi+fi- 


-1.3 ±2.8 ±1.1 


0.65 


16 


4.3 


D+ -> K+e+^i- 


-4.3 ±1.8 ±0.6 


1.44 


4.3 


1.2 


D+ -^ K+fi+e- 


3.2 ±3.8 ±1.2 


1.74 


9.9 


2.8 


D+ -> K+e+e- 


-5.7 ±5.8 ±2.0 


3.20 


1.6 


3.7 


D+ -^ K+ti+fi- 


4.8 ±5.9 ±1.2 


0.85 


9.1 


21 


D+ -^ K+e+^jL- 


9.1 ±6.0 ±2.8 


1.74 


5.7 


14 


Dt -^ K+^ji+e- 


3.4 ±6.4 ±3.5 


2.08 


4.2 


9.7 


At -> pe+e" 


4.0 ±6.5 ±2.8 


5.52 


0.8 


5.5 


At -^ PM"^M" 


11.1 ±5.0 ±2.5 


0.86 


6.4 


44 


At -> pe+^- 


-0.7 ±2.9 ±0.9 


1.10 


1.6 


9.9 


At — > pfj.'^e' 


6.2 ±4.6 ±1.8 


1.37 


2.9 


19 


D+ ->7r-e+e+ 


4.7 ±4.7 ±0.5 


3.16 


6.8 


1.9 


D+ ->7r-^+/i+ 


-3.1 ±1.2 ±0.5 


0.70 


7.5 


2.0 


D+ -^7v-fi+e+ 


-5.1 ±4.2 ±2.0 


1.72 


7.4 


2.0 


Dt -^Tv-e+e+ 


-5.7 ±14. ±3.4 


6.84 


1.8 


4.1 


Dt -^ 7r-^+/i+ 


0.6 ±5.1 ±2.7 


1.05 


6.2 


14 


Dt -^ Tr'fi+e+ 


-0.2 ±7.9 ±0.6 


2.23 


3.6 


8.4 


D+ -> K-e+e+ 


-2.8 ±2.4 ±0.2 


2.67 


3.1 


0.9 


D+ ^ K-fi+^i+ 


7.2 ±5.4 ±1.6 


0.80 


37 


10 


D+ ^ K-fi+e+ 


-11.6 ±4.0 ±3.1 


1.52 


6.8 


1.9 


Dt -^ K-e+e+ 


2.3 ±7.9 ±3.3 


4.10 


2.1 


5.2 


Dt -^ K-fi+fi+ 


-2.3 ±5.0 ±2.8 


0.98 


5.3 


13 


Dt -^ K-fi+e+ 


-14.0 ±8.4 ±2.0 


2.26 


2.4 


6.1 


At -> pe+e+ 


-1.5 ±4.2 ±1.5 


5.14 


0.4 


2.7 


At -^ m+iJ.+ 


-0.0 ±2.1 ±0.6 


0.94 


1.4 


9.4 


At -^ pfj.'^e'^ 


10.1 ±5.8 ±3.5 


2.50 


2.3 


16 



VIII. CONCLUSIONS 



Searches for the decay modes D 



Dt, -^ K^l^&^+ and A+ -^ 

is) c 

formed using 384 ftj- of e+e 



is) 



7r±^=F^(')+, 



t)£T^(')H 



pt^V ' have been per- 

annihilations collected 
with the BaBAR detector. No signals are observed and 
we report upper limits on 35 different branching ratios 
between 0.4 x 10"'' and 37 x lO"'^ at 90% CL. This cor- 
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FIG. 4: Invariant-mass distributions for D"*" ^- tt'^^^^''^^ (top) and Df -^ 7r+^+^'''^ (bottom) candidates. The solid lines are 
the results of the fits. The background components for the dimuon modes in which muon candidates arise from misidentified 
hadrons are shown as dashed curves. 
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FIG. 5: Invariant-mass distributions for D"*" — ^ K'^£^v''~ (top) and Df — > K'^i^v''" (bottom) candidates. The solid lines are 
the results of the fits. The background components for the dimuon modes in which muon candidates arise from misidentified 
hadrons are shown as dashed curves. 



responds to limits on the branching fractions between 
1 X 10~^ and 44 x 10~^. These limits are calculated un- 
der the assumption of three-body phase-space decays; the 
efficiency varies by up to 25% as a function of dilepton 
invariant mass. For 32 of the 35 decay modes studied, 
the limits are an improvement over the existing measure- 



ments and therefore provide more stringent constraints 
on physics beyond the SM. 
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FIG. 6: Invariant-mass distributions for I?"*" ^- tt"^"*"^ '■''"'" (top) and Df -^ iv~t^t'''^ (bottom) candidates. The solid lines are 
the results of the fits. The background component for the dimuon modes in which muon candidates arise from misidentified 
hadrons is shown as a dashed curve. 
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FIG. 7: Invariant-mass distributions for D+ ^ A'"£+£<''+ (top) and D+ ^ K' 1+ &^+ (bottom) candidates. The solid lines 



are the results of the fits. The background components for the dimuon modes and DJ, 
from misidentified hadrons are shown as dashed curves. 



K fj, e in which candidates arise 
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FIG. 9: Invariant-mass distributions for A^ — >■ pi'^ v-''^ candidates. The solid lines are the results of the fits. 
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FIG. 



10: Invariant mass distribution for (left) Dt. 



e+e- and (right) D+j -> tt"* 



candidates. Only 



events with dilepton invariant masses in the (f> signal region 
defined in Sect. IIV Al are plotted. The solid lines are the re- 
sults of the fits. 



TABLE V: Signal yields for the fits to the Dt. -^ tt^(J>i+i- candidates. For the yields, the first uncertainty on the yield is 
statistical and the second is systematic. The last column shows the expected number of signal events, where the uncertainty is 
due to the systematic uncertainty assigned to the efficiency. 



Decay mode 



Yield (events) 



Efliciency (%) 



Expected yield (events) 



D+ 
D+ 

Dt 
D+ 



"e+e- 
t>e+e- 



21.8 ±5.8 ±1.5 

7.5 ±3.4 ±1.4 

63 ± 10 ± 3 

12.7 ±4.3 ±2.6 



5.65 
1.11 
6.46 
1.07 



22.2 ±1.1 
4.5 ±0.4 
79 ±3 
13.1 ±1.2 



